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^ ; ABSTRACT 

^< ' We study the structure and evolution of "quasistars," accreting black holes embedded 

"^h \ within massive hydrostatic gaseous envelopes. These configurations may model the 

\0 ■ early growth of supermassive black hole seeds. The accretion rate onto the black hole 

adjusts so that the luminosity carried by the convective envelope equals the Eddington 
limit for the total mass, M* + A/bh ~ M*- This greatly exceeds the Eddington limit 
for the black hole mass alone, leading to rapid growth of the black hole. We use 
analytic models and numerical stellar structure calculations to study the structure and 
evolution of quasistars. We show that the photospheric temperature of the envelope 
scales as T p h oc M BH MJ , and decreases with time while the black hole mass 
increases. Once T p h < 10 4 K, the photospheric opacity drops precipitously and T ph 
hits a limiting value, analogous to the Hayashi track for red giants and protostars, 
below which no hydrostatic solution for the convective envelope exists. For metal-free 
(Population III) opacities this limiting temperature is approximately 4000 K. After 
a quasistar reaches this limiting temperature, it is rapidly dispersed by radiation 
pressure. We find that black hole seeds with masses between 10 3 M© and 10 4 M© 
could form via this mechanism in less than a few Myr. 
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1 INTRODUCTION 

The formation mechanism for supermassive black holes 
in galactic nuclei remains unknown. V ariations on most 
of the formation channels identified by iBegelman fc Reesl 
1)19781 ) which include instabilities in clusters of 

stars or stellar remnants, and the collapse of su- 
permassive stars or massive discs are still 
der consideration 



today |Umemura. Loeb fc Turnen 
Rasi j j2006l : IShibata fc Shapirol 



1993 



Freitag, Gurkan fc Rasiol 120061 : IShibata fc Shapirol 12002 
Lodato & Nataraianj [20061 ). A more recent idea holds 
that supermassive black holes result from sustained ac- 
cretion onto, or mergers of, the remnants of Popula- 
tion III (i.e. metal-free) stars (|Volonteri. Haardt fc Madaul 
I2003T ). at least some of which seem likely to be 
massive, short-lived progenitors of stellar mass black 



holes dCarr. Bond fc Arnettlll984 
19991: 



, Abel, Br yan, fc Norma: 
Tumlinson. Venkatesan fc ShulJ[20ol ) 



Bro mm. Coppi fc Larson! 
20021: iHeger et al.l |2003I 



No current observation directly constrains the different 
classes of model for black hole formation. However, there 
are clues. The existence of massive black holes in quasars at 
redshifts z > 6 strongly suggests that these black holes, at 



least, started forming at redshifts high enough (z ~ 20) to 
predate extensive metal enrichment by the first stars. Moti- 
vated by this, we develop in this pa per a model proposed by 
IBegelman. Volonteri fc Reesl (|2006T ) (BVR), in which super- 
massive black holes form, not from Pop III stars themselves, 
but rather from the evolution of a new class of Pop III ob- 
jects that might form in metal-free haloes too massive to 
yield individual stars. In outline, the BVR model envisages 
a three-stage process for black hole formation. First, gas in 
metal-free haloes with a virial temperature above T ~ 10 4 K 
flows toward the center of the potential as a result of gravi- 
tational instabilities, forming a massive, pressure-supported 
central object. Nuclear reactions may start, but the very 
high infall rate continues to compress and heat the core, 
precluding formation of an ordinary star. Eventually, when 
the core temperature attains T~5x 10 s K, neutrino losses 
result in a catastrophic collapse of the core to a black hole. 
We dub the resulting structure — comprising an initially 
low mass black hole embedded within a massive, radiation- 
pressure supported envelope — a quasistar. Initially, the 
black hole is much less massive than the envelope. Over time, 
the black hole grows at the expense of the envelope, until 
finally the growing luminosity succeeds in unbinding the en- 
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Figure 1. Schematic illustration of the quasistar structure that 
we consider in this paper. A seed black hole of mass Mgu ac- 
cretes gas from a massive, radiation pressure-supported envelope 
at a rate set by the conditions outside the Bondi radius. The 
luminosity liberated by the accretion process is transported con- 
vectively in the inner regions of the envelope, with a transition to 
a radiative zone once convection becomes inefficient. In this paper 
we consider isolated, spherically symmetric models of quasistars, 
but the more physical situation would also include partial rota- 
tional support leading to flattening of the quasistar, and ongoing 
disc accretion at a fraction of a Solar mass per year. 



velope and the seed black hole is unveiled. The key feature of 
this scenario is that while the black hole is embedded within 
the envelope, its growth is limited by the Eddington limit 
for the whole quasistar, rather than that appropriate for the 
black hole mass itself. Very rapid growth can then occur 
at early times, when the envelope mass greatly exceeds the 
black hole mass. 

The structure of quasistars, illustrated in Fig. [l] re- 
sembles that of more familiar objects. The outer regions 
have qualitative similarities to Thorne-Zytkow (1977) ob- 
jects, although in our case the luminosity derives entirely 
from black hole accretion. The accretion p hysics is related 
to the "hyper critical" regime of accretion (IBegelmanl 1 19791 ; 
iBlondinl Il986f ) , which may be observed in X-r ay binaries 
such as SS433 (|Begelman, King fc Pringiel 12006ft . However, 
these are very rough analogues. The large difference between 
Pop III opacity and that of present-day, metal-enriched gas, 
together with the dominant role of radiation pressure and 
the presence of ongoing accretion in quasistars, result in a 
unique quasistar structure. 

In this paper, we study the structure and evolution of 
quasistars. Using Pop III opacities, we construct simplified 
analytic models which we compare with numerical integra- 
tions of the stellar structure equations. Our model quasistars 
are assumed to be spherically symmetric and in hydrostatic 
equilibrium. Real quasistars, if they exist, may not obey 
these strictures, and thus our results cannot be definitive. 
Rather, our calculations are intended to provide a first esti- 
mate of the maximum black hole mass that this mechanism 
can yield, together with a determination of the photospheric 
temperature and luminosity of quasistars that can be used 
to assess their detectability with next-generation observato- 
ries such as the James Webb Space Telescope. 



2 ANALYTIC CONSIDERATIONS 

We consider a black hole of mass Mbh embedded within an 
envelope of mass M* 3> A/bh- The luminosity of the quasis- 
tar is generated exclusively by black hole accretion, whose 
rate depends upon the central conditions in the envelope. 
To proceed analytically, we first develop simple scaling re- 



lations for the envelope structure and black hole accretion 
rate. 



2.1 Envelope Structure 

Quasistar envelopes with masses > 10 3 Mq are supported 
primarily by radiation pressure. Since the luminosity car- 
ried by the envelope must equal the Eddington limit for the 
total mass, the flux carried by radiative diffusion in the en- 
velope's interior is only a fraction M(< r)/M» of the total, 
where M(< r) is the mass enclosed within r. Thus, quasis- 
tar envelopes are strongly convective, and their structures 
resemble n — 3 (7 = 4/3) polytropes. The most accurate 
approa ch would be to model the envelope as a "loaded poly- 
trope" IjHuntlev fc Saslawill975r ) , with the black hole treated 
as a central point mass, but for M* 3> Mbh the standard 
Lane-Emden solutions suffice. Defining m* = M*/Mq, the 
ratio of gas to radiation pressure is uniform throughout the 
convective zone. We find 

Ba = 7.1m»- 1/2 , 

Pr 



(1) 



where we have assumed a mean mass per particle p (a 0.6m p 
valid for fully ionized regions. 



2.2 Energy Source 

The central regions of n = 3 polytropes have approximately 
uniform densities (p c ), temperatures (T c ), and pressures 
(p c — aTc/3). The boundary conditions for black hole accre- 
tion are therefore similar to those of Bondi (1952) accretion, 
or more specifically the gen eralization for optical thick flows 
given by iFlammangl |l982j). The adiabatic accretion rate is 



Mb 



4tt (GM B h) 2 /Qc 
V2 d 



(2) 



where c c = (Ap c /'ip c ) 1 ' 2 is the adiabatic sound speed. 
Bondi's solution assumes that all gravitational binding en- 
ergy liberated during accretion is advected into the hole, but 
this is unrealistic in the presence of even a small amount of 
rotation. Provided that the specific angular momentum at 
the Bondi radius Zfl ow exceeds the specific angular momen- 
tum Z ms of the marginally stable circular orbit, we expect 
the flow to be rotationally supported. In this case a geo- 
metrically thick accretion disk will form around the hole, 
within which angular momentum transport is required in 
order for accretion to occur. Although the efficiency of such 
a disk (and how the energy output is partitioned between 
radiation and mechanical work) is hard to calculate, it is 
reasonable to assume that it will depend primarily on the 
depth of the black hole potential. We write the luminosity 
as Lbh = eA/bhc 2 , where e ~ O(0. 1) is the efficiency of 
energy output and A/bh is the actual accretion rate. For 
simplicity, we take e to be a constant. In the absence of 
an efficient exhaust such as a jet or evacuated funnel, this 

1 We note that e is expected to vary by factors of the order of 
unity with the black hole spin. In our models, however, we are 
envisaging a black hole that grows in mass by several orders of 
magnitude. If the accretion flow retains a similar character during 
this growth, we expect that the spin of the hole will rapidly attain 
a limiting value, after which further changes will be small. 
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energy must be carried beyond the Bondi radius convec- 
tively, by the accreting gas itself. Given that the convective 
flux density may not exceed ~ PcC c , we conclude that the 
accretion rate will be reduced below the Bondi value by a 
factor of order e _1 (c c /c) 2 <C 1 (Gruzinov 1998; Blandford 
& Begelman 1999; Narayan, Igumenshchev & Abramowicz 
2000; Quataert & Gruzinov 2000). Unless e < 0.1 - for ex- 
ample in the case where Zfl ow < l ms - accretion is suppressed 
by a large factor, of the order of 10 4 for the models we con- 
sider. Combining these results we adopt the expression 



47rG 2 aM BH p 3/2 p- 1/2 



(3) 



for the accretion luminosity. The parameter a < 1 accounts 
for energy sinks within the Bondi radius: inefficient convec- 
tion, presence of outflows, etc.; as well as any inefficiency 
of angular momentum transport. Small values of a imply a 
reduced energy supply to the quasistar. The standard result 
for ann = 3 polytrope (e.g. Hoyle & Fowler 1963), 



(4) 



p c = 1.3 x 10~ 4 m7 1/2 T 6 3 g cm" 3 , 

where Ts = T c /10 6 K, coupled with the equation of state for 
radiation in LTE, allows us to express Lbh in terms of the 
central temperature: 



a ., 1n 42 2 -3/4^5/2 

'.6 x 10 am Bn m t ' i 6 erg s 



(5) 



where ttibh = Mbh/M^ 



2.3 Scaling Laws 

Suppose that the convective zone encompasses nearly the 
entire mass and radius of the envelope, allowing us to treat 
M* and R* as constant within the radiative layer. We will 
check the validity of this assumption in § 2.4.2. We expect 
Lbh to be very close to the Eddington limit at the transition 
between the convective zone and the outer layer where radia- 
tive diffusion carries all the flux, and set Le = 4 7rGM,c/ /ttr, 
where «t r is the opacity at the transition radius. As we will 
see below, Kti is close to the electron scattering opacity, 
k ce = 0.35 cm 2 g _1 , hence k = Ktr/^cs ~ 0(1). We therefore 
write 



L BH = 1.4 x 10 38 '^k' 1 m* erg s" 1 , 



(6) 



where £ tI w 1 is the Eddington factor at the transition. 
Equating the two expressions for Lbh, we obtain 



T c = 1.4 x 10 X 



4 ,,2/5 ~ -2/5 -2/5 -4/5 7/10 



K. 



(7) 



We are interested in systems with m, ~ O(1007tibh), so T c 
typically lies in the range 10 5 — 10 6 K and electron scattering 
opacity dominates in the quasistar interior. At these temper- 
atures, the rates of energetically significant thermonuclear 
reactions are negligible, and can be safely ignorecQ We note 
that this differs from the case of envelope accretion onto 
neutron stars, where the presence of a hard surface results 
in high temperatures and potentially significant nuclear re- 
actions within the hydrostatic region of the flow. Although 

2 In the most massive quasistars, the central temperature may be 
high enough (a few million K) to initiate lithium burning. This is 
energetically negligible, and although the presence or absence of 
lithium does affect the opacity, the effect is small for the photo- 
spheric temperatures and densities of interest here. 



we do not model the non-hydrostatic region of quasistars 
in any detail, in our case high temperatures are attained in 
the immediate vicinity of the black hole. However, even in 
this region the neglect of nuclear reactions is justified, first 
because black hole accretion is energetically much more ef- 
ficient than fusion, and second because ongoing accretion 
limits the timescale over which inflowing gas is exposed to 
high T. 

The polytropic relations also allow us to estimate the 
radius of the quasistar, 

R* = 5.8 x 10 12 my 2 r 6 _1 cm, (8) 

l|HovIe fc Fowlerlll963r ) and therefore its photospheric (effec- 
tive) temperature, 

T ph = 8.5 x 10 3 £ t 1 / 4 £- 1/4 T 6 1/2 K. (9) 

Inserting our estimate for T c , we obtain 

R. = 4.3 x 10 14 ^- 2/5 S 2/5 a 2/5 m 4/ H 5 m; 1/5 cm (10) 

and 



T ph = 1.0 x 10°C k 



3^9/20^-9/20^-1/5^-2/5^7/20 R 



m* 



(11) 



Thus, quasistars should have radii of order 10 2 — 10 3 AU 
and temperatures of a few thousand degrees. 



2.4 Radiative Layer and Photosphere 

In the outer layers of the quasistar, convection is unable to 
transport the total luminosity and radiative diffusion takes 
over as the dominant energy transport mechanism. To esti- 
mate the transition temperature between the convective and 
radiative zones, we note that the maximum flux that can be 



transported convectively is 



L 



con, max 



(3pc s 



(12) 



where c s = (p/p) 1 ' 2 is the local (isothermal) sound speed 
and /3 < 1 is an efficiency factor. Equating this to the pho- 
tospheric flux acT 4 h /4 (again assuming a narrow radiative 
zone), we obtain 



T 4 



3 C 4 



4/3 c s 



pi, • 



(13) 



This result is equivalent to the condition that radiation be 
trapped in the convective cells as they rise. Since j3 < 1 and 
c/c a 2> 1, we deduce that T tr 3> Tph- Expressing p and p in 
terms of T tI (using the polytropic relations and the equation 
of state), we obtain 



T tr = 31/3 



-2/9^-1/18718/9 



T. 



P h 



(14) 



Since we will later show that T p h cannot drop below a few 
thousand degrees, we conclude that T tr is well above 10 4 
K. The import of this is that at the densities and temper- 
atures likely to apply near the transition (for m* > 10 4 ), 
bound-free opacity can be important but will not elevate 
the Rosseland mean opacity above Thomson scattering by 
a large factor. 

In order to determine the structure of the radiative 
layer, we need to know the behavior of the opacity. For the 
purpose of our analytic estimates, we adopt a simple phe- 
nomenological form for the opacity, based on the Pop III 
opacity tables of Mayer & Duschl (2005, hereafter MD05). 
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We note that, for temperatures below 10 4 K and densities 
< 10~ 9 g cm -3 , the Rosseland mean opacity depends much 
more sensitively on temperature than on density. This can 
be seen clearly in Fig. 4 of MD05. We therefore write 



«(T) = 



Ml 



1 + (T/T ) 



(15) 



This expression does not capture the possible contribution 
of bound-free opacity above 10 4 K, but it does mimic the ex- 
tremely steep decline in opacity toward lower temperatures. 
A crude fit gives s « 13, To « 8000 K, and kq « k cb . 

The advantage of using an opacity that is solely a func- 
tion of temperature is that one can combine the equation of 
hydrostatic equilibrium with the radiation diffusion equation 
and integrate to obtain the total (radiation + gas) pressure 
as a function of temperature (see, e.g., Cox 1968, chapter 
20, or any stellar structure textbook for a discussion of this 
technique). Defining the Eddington factor associated with 
opacity kq as 



Lbhko 
47tGM*c' 



we obtain 

4a 

Wo 



p(T) = 



T 1 



Ti4 
ph 



+ 



To(TX s -T 4 - 



+ Pph 



(16) 



(17) 



where p p h is the total pressure at the photosphere. We 
use a standard model for the atmospheric structure, in 
which the radiation pressure depends on optical depth as 
Prij) = (aT ph /6)(l + 3T/2) and the temperature is constant 
between the true (r = 0) surface and the photosphere at 
r = 2/3 (see, e.g., Mihalas 1978 or most stellar structure 
textbooks). Integrating the equation of hydrostatic equilib- 
rium and noting that the surface gravity can be written as 
g = KoaTp h /(4:£o), we have 



Pph = r 



3 K(Tph) 



+ ■ 



aT, 



P h 



aT, 



P h 



(i 



6/„ 



Jph 



(18) 



P S (T) 



The gas pressure is then given by 
aT 4 



+ 



4:0, 

34 


/ T 4 


-t) 


TS 

s — 


/s + 4 4 _ s 
4 v 8 ph 


_ T 4-sY 



(19) 



for temperatures T p h ^ T ^ T tI . 

The ratio of gas pressure to radiation pressure must be 
continuous across the boundary between the convective zone 
and the radiative layer. The matching condition is, therefore, 



3p 9 (T tr 
aT 4 



= 7.1m; 



1/2 



(20) 



To simplify matters, we note that, since Ttr 3> T pri and 
s ^> 4, we can neglect the terms in p g (T) proportional to 



'ph 

obtain 



T„ h and T s . Solving the matching condition for £q, we 



i+ 



<?0 = 



s+4 
2(s-4) 



1 + 7.1ml 



1/2 



(21) 



2.4-1 Opacity Crisis 

A necessary condition for a model to be physically realistic 
is that the flux not exceed the Eddington limit at the tran- 
sition radius, £ tI < 1. At the transition radius, radiation is 
effectively trapped, and any super-Eddington flux would be 
efficiently converted into bulk kinetic energy. Hence, if this 
condition is violated we expect the outer layers to be blown 
off by radiation pressure, ultimately dispersing the quasis- 
tar. For our approximate opacity, which is a monotonically 
increasing function of T, this translates to £o < l + (To/Tt Y ) s . 
Since the convective envelope is dominated by radiation 
pressure, we can write this condition as 



s + 4 
2(s - 4) 



(!±Y (Ik 



< 7.1m,: 



- 1/2 + (g) . ,22, 



Noting that the second term on the right-hand side of this 
equation is negligible, we substitute for T tr using equa- 
tion (|14[) to obtain a lower limit on T p h, 



t;- 4/9 > 7.7 x 10~ 8 



(£)' 



8/9 13/18 T » 

m« Jq ■ 



(23) 



For s>4, this lower limit is very insensitive to the convec- 
tive efficiency j3 and the envelope mass M*, and is roughly 
proportional to Tq. For example, adopting s — 13 and nor- 
malizing /3, m* and To to 0.1, 10 , and 8000 K, respectively, 
we have 



,8/113 13/226^,117/113 
' r "*4 1 0, 8000 



Tph > 4500 fjj_ t m^ A ' T n a ' nnn K. 



(24) 



The corresponding lower limit on T tr shows similar behavior: 

T v re nnn 0-18/113^-1/226^104/113 v /„.v 

T tr > 55,000 p_ 1 m, 4 ' T ^ 000 K. (25) 

The first of these lower limits is the most important 
result of the analytic part of this paper. It represents a floor 
to the photospheric temperature of quasistars, analogous to 
the "Hayashi track" (Hayahi & Hoshi 1961; Hayashi 1961), 
which limits the temperatures of red giants and convective 
protostars. Our analysis differs from Hayashi's model in that 
our convective envelopes are radiation pressure-dominated 
and therefore resemble n — 3 polytropes, rather than the 
n = 3/2 polytropes appropriate to gas pressure-dominated 
convection. 



2.4-2 Validity of Approximations 

The analytic model described above assumes that the 
growth of the black hole within the envelope can be 
described by a sequence of hydrostatic solutions. This 
quasi-static approximation is justified by the ordering of 
timescales. Adopting representative values of the mass, ra- 
dius, and radiative efficiency (m» = 10 5 , jtibh = 10 3 , 
R, — 10 15 cm, e = 0.1) the timescale on which hydrostatic 
equilibrium is established, 



£dyn 



GM„ 



R* 



10 B s 



(26) 



is much shorter than the timescale on which the black hole 
grows 

Mbh 



Mbh 



10 13 s. 



(27) 



Similarly the timescale for attaining radiative equilibrium, 
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tr 



rAR, 



10 7 S, 



(28) 



where r is the optical depth through the radiative layer of 
width A_R», is enormously shorter than any of the evolu- 
tionary timescales of the system. 

In addition to the quasistatic assumption, we have taken 
the geometric thickness and mass of the radiative layer to be 
negligible compared to R* and M, , respectively. How good 
are these assumptions? 

First, consider the mass of the radiative layer. If the 
layer is geometrically thin and the gravity g is approximately 
constant across it, then we have 



AM. _ 4TTRJaT^ 
M* ~ 3GM„ 2 



(29) 



from the equation of hydrostatic equilibrium. Using equa- 
tions p0|) . (Ilip . and (|14p to eliminate _R„, T tI , and to,bh in 
favor of T p h, we obtain (for £ tr = ft = 1) 



AM. 
M. 



0.2/3. 



-8/9 -2/9 



V4500 K7 



-40/9 



(30) 



where /3-i = /3/0.1. Equation (|30|) indicates that the as- 
sumption of constant mass in the radiative zone is only 
marginally self-consistent when T p h is close to its minimum 
value. The approximation improves at larger T p h. 

The geometric thickness of the radiative layer is domi- 
nated by the region close to the transition temperature. The 
opacity is therefore close to kq ~ k cs and we may approxi- 
mate 



A_R» 
R* 



167vR,acT^ 

OUcsptr 



0.7/3. 



-2/9 



-1/M / 



tph 



V4500 K 



-10/9 



(31) 



Evidently, assuming geometrical thinness in the radiative 
layer is an even poorer approximation than assuming con- 
stant enclosed mass. In fact, our numerical models show 
that even the scaling in eq. (|31[) - obtained assuming 
AR*/R* < 1 — is not valid. The radiative layer thickness 
generally decreases with T p h. 



Fi 



4 r 4 I p g 
3 3 6 p 



(33) 



If the pressure- weighted integral over the star of 57 < Sy CI it, 
instability is possible. 

The region of the quasistar interior to the Bondi radius 
is not in hydrostatic equilibrium, so it is meaningless to ap- 
ply the above condition there. For an estimate, we evaluate 
the above expressions at the Bondi radius using the expres- 
sions for the interior structure of the quasistar. We find that 
the change in the stability boundary is 



4 x 10- 10 m! /5 m-^ V 



2/5 



*7dt = (^) 

while the actual 7 exceeds 4/3 by an amount 



57 ~ 1.2m : 



-1/2 



(34) 



(35) 



We conclude that instability is possible if m* is large enough. 
For a — 0.01 and rriBH = 10 3 , for example, instability is 
possible for m* > 3 x 10 6 . 

We caution that this analysis is far too simple to de- 
scribe the actual situation. Globally, a robust expectation 
is that quasistars will be rapidly rotating as a consequence 
of their formation from rotationally supported gas (BVR). 
Unfortunately, this does not help us in estimating the ro- 
tation rate at the Bondi radius, as the latter depends upon 
the very uncertain role of convection in redistributing an- 
gular momentum. However, any significant rotation would 
help stabilize the structure against dynamical instability. 



3 NUMERICAL MODELS 

Although our analytic models ought to yield a reasonably ac- 
curate picture of the generic features of quasistar envelopes, 
they are only marginally self-consistent for photospheric 
temperatures approaching the floor value. Moreover, real 
Pop III opacities, which are functions of density as well as 
of temperature, have considerably more structure than the 
analytic fit employed in the analytic models, and require a 
numerical treatment. 



2.5 Stability 

The interior conditions of quasistars are hot enough that the 
opacity is dominated by electron scattering, yet too cool for 
nuclear reactions to occur. Accordingly, we do not expect 
any of the stellar instabilities that depend upon complex 
opacities or nuclear reaction rates to afflict the interior of 
quasistars (we note later the possibility of an instability near 
the surface due to the presence of locally super-Eddington 
zones). A more serious concern - given that quasistars are 
highly radiation pressure dominated - is dynamically insta- 
bility. In a non-rotating model, relativistic effects raise the 
critical 7, below which instability occurs, to 



7crit 



- + 5-Jcrit 



(32) 



where S^u ~ GM/Rc 2 (e.g. Shapiro & Teukolsky 1983). 
This must be compared to the actual 7 in the radiation 
dominated interior of the quasistar. For p g /p <C 1, the first 
adiabatic exponent is given by 



3.1 Equations 

We model quasistars as static, spherically symmetric objects 
in thermal equilibrium. The equation of hydrostatic equilib- 



dp G M{r) 

T = 2 — P' 

dr r A 

where the mass enclosed within radius r is given by 

dM 
dr 



■ 4nr p. 



(36) 



(37) 



At the Bondi radius, which for our numerical integrations 
we define as 



Rl: 



2p c ' 



(38) 



we set M(i?Bondi) = Mbh- The gas mass within the Bondi 
radius is negligibly small compared to the black hole mass, so 
ignoring the non-hydrostatic region in our structure calcu- 
lations should be an excellent approximation. The equation 
of state is a mix of gas and radiation pressure, 



© 0000 RAS, MNRAS 000, 000-000 



6 Begelman, Rossi & Armitage 

p = Pa +Pr = ^ + \aT\ (39) 

/i 3 

where k is the Boltzmann constant and p the mean molecu- 
lar weight. We set p = 0.6m p , appropriate for a fully ionized 
gas, and ignore the variation of p due to partial ionization 
near the quasistar photosphere. This is a very good approx- 
imation, since the neutral hydrogen fraction generally does 
not exceed a few per cent. 

The temperature gradient is determined by the relative 
amount of energy carried by convection (F cou ) and radiation 
(-frad) at each radius. Since quasistars are powered solely 
by black hole accretion, the luminosity Lbh = 47rr 2 (.F con + 
-Frad) is constant with radius. To determine the boundary of 
the convective zone, we use the Schwarzschild criterion for 
convective stabilitjQ 

dT ra a dT a d 



>0, 



(40) 



dr dr 

where the radiative gradient dT ra< j/dr is the gradient of tem- 
perature if all the energy is transported by radiative diffu- 
sion, 



rfTrad 3_ Lbh pn 

dr Aac 4nr 2 T 3 ' 

and the adiabatic gradient, 

dTkd = / r 2 -l \ Tdp 
dr \ Y-2 J p dr' 



(41) 



(42) 



describes the temperature variation in a convective element 
as it moves adiabatically (no radiative losses) through the 
surrounding layers. The adiabatic exponent F 2 is given by 



r 2 = 



32 - 24/3 - 3f3 2 
24 - 18/3 - 3f3 2 



(43) 



where /3 — p g /p is the fraction of the total pressure con- 
tributed by gas pressure. We then set the actual temperature 
gradient to be, 

dTrad 1 

i, (44) 



dT 
dr 



dT ri 



[min(- 



dr 



+ 



dr 



where x — F/F cont max . Typically, we find that the pressure 
scale height (which in a mixing length theory for convection 
is comparable to the distance convective elements travel) is 
a small fraction of the radius across the outer «90% of the 
quasistar by radius. In the inner regions, where the pressure 
profile is quite flat, x is very small, and hence it is reason- 
able to assume that convection (even if it behaves differently 
from a mixing length theory) also succeeds in establishing 
the adiabatic gradient there. As in the analytic models, we 
note that the maximum flux that convection can transport 
is limited by the condition that convective motions cannot 
be supersonic. Accordingly, we set 



F« 



con, max 



: f3c s p r 



(45) 



with the efficiency parameter )3 = 0.1. 

Equation (|44p describes the limiting behaviors of the 
temperature gradient in the different regions of a quasistar. 

3 Note that, in using the Schwarzschild rather than the Ledoux 
criterion for convective stability, we again ignore effects due to 
changes in the mean molecular weight with radius. As explained 
above, we believe that this is a good approximation. 



The adiabatic gradient applies if the region is dynamically 
unstable and the flux is less than the maximum convective 
flux; otherwise, the radiative gradient applies. These limit- 
ing cases accurately model the regions of high convective 
efficiency, in the denser and more opaque interior layers, 
and of high radiative efficiency, in the generally stable outer 
layer. In the transition region, however, our use of a flux lim- 
iter (the 1 + a; 10 term in eq. [44]) is only an approximation. 
A more accurate treatment would make use of mixing length 
theory. However, we have checked that different descriptions 
of the transition region do not affect appreciably the overall 
quasistar structure. 



3.2 Integration Method 



Equations (J36J) , J37J, (J39J) , and (Q4J), together with the defi- 
nitions of the radiative and adiabatic temperature gradients, 
form a closed set of equations for the 4 unknowns (p, T, p 
and mass) , given a luminosity Lbh and stellar radius R, . To 
solve this system, we fix the photospheric temperature T p h, 
the black holes mass Mbh, and the parameter a, and guess 
the photospheric radius R, and the quasistar mass M*. The 
luminosity is then L 8 g = AixR^aT^, where a — ac/4. We 
determine the photospheric pressure and density using the 
first equality in eq. (|18[1 , together with the equation of state 
(eq. [39]). We then integrate inward from the known pho- 
tospheric conditions, adjusting our initial guesses until we 
match the two conditions 



Left = Lbh (10 X ifeondi) 
M»(iWdi) = 0. 



(46) 
(47) 



We calculate Lbh (eq. [3]) with p and p evaluated at 
10 x -RBondi, well beyond the inner region affected by the 
black hole gravity. For our numerical examples a is assumed 
to equal 0.1, if not otherwise stated. We discuss the depen- 
dence of our results on a in § 14.31 We explore two models for 
the opacity: the "toy" opacity given by eq. (|15[1 . which we 
used previously in the analytic models, and the full numeri- 
cal opacity tabulated by MD05, which we refer to as "Pop III 
opacity" . Fig. [2]shows a comparison of these opacities for dif- 
ferent densities of interest in the outermost radiative layers. 



3.3 "Toy" Opacity 

As a test, both of our analytic models and of our numer- 
ical scheme, we first compute structure models using the 
"toy" opacity (eq. I15[) . The toy opacity ignores density de- 
pendence, and in particular the contribution from bound- 
free and free-free absorption at T > 8 x 10 3 K. As can be 
seen in Fig. [2] this is a reasonable approximation at low 
density, but is very poor at higher density. 

As expected, we find that the analytic scalings are most 
reliable for high photospheric temperatures, and less reliable 
close to the temperature floor. We confirm that for a given 
(a, Mbh) there is indeed a minimum photospheric temper- 
ature, T m i n , above which the luminosity at the transition 
radius is sub-Eddington. Below T m i n , the whole radiative 
layer experiences a radiative force greater than the gravi- 
tational force and the quasistar is bound to evaporate. In 
Fig. [3] the short-dashed line shows T m i n as a function of 
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Figure 2. Opacity versus temperature. The Pop III opacity from 
MD05 is plotted as a function of temperature for three differ- 
ent values of matter density as shown in the legend. The range 
10~ 13 g cm -3 $C p ^ 10~ 8 g cm -3 spans the densities to be 
found in the outermost layers of quasistars. The solid line shows 
the analytic opacity given by equation H15H • The figure shows how 
the Pop III opacity increases over the analytic fit as the density 
increases. The effect is especially evident around T = 10 4 K, 
where there is a bound-free peak due to hydrogen ionization. 



Mbh for numerical models computed using the toy opac- 
ity. The analytic estimate, obtained by combining equations 
Q24p and pip , is shown as the long-dashed line. Numerically, 
we obtain a slightly higher T m i n . The discrepancy is higher 
for higher black hole masses. We note that for models com- 
puted using the toy opacity, the critical temperature below 
which the entire radiative layer becomes super-Eddington 
(our definition of Tmin) is close to the temperature at which 
any point in the radiative layer becomes super-Eddington. 
This simple behavior, which occurs because the toy opacity 
is constant in most of the radiative layer before dropping off 
steeply and monotonically near the photosphere, does not 
carry over to models computed using full Pop III opacity. 



3.4 Pop III Quasistars 

To extend these results, we compute numerical quasistar 
structures with realistic opacities. We assume that quasis- 
tars have a primordial (metal-free) nuclear composition, and 
use the opacity table calculated by MD05. This table covers 
the density range 1CP 16 < p [gem -3 ] < 1CP 2 for tempera- 
tures 63 < T [K] < 4x 10 4 . To compute the interior structure 
of the quasistar — where temperatures attain values much 
higher than 4 x 10 4 K — we analytically extend the opacity 
table assuming that the excess opacity above the electron 
scattering value has a Kramers form, i.e., that, 



K(p,T) = C(p) P T- 



+ K e 



(48) 



We fix the function C(p) to match smoothly onto the tabu- 
lated opacity at the highest temperature point. As expected, 
Thomson scattering is the only significant source of opacity 
deep within the interior. 

As is evident from our Fig. [2] or from Fig. 4a of 
MD05, the opacity deviates substantially from our toy model 




1000 



Figure 3. Minimum photospheric temperature T m ; n versus black 
hole mass in unit of solar masses mBHi f° r a = 0.1. Numerical 
models with the "toy" opacity (short-dashed line) are compared 
to results computed using the Pop III opacity (solid line). The 
analytic estimate, obtained by combining equations (1241 and llllll . 
is shown as the long-dashed line. Static quasistar models do not 
exist in the lower shaded region. 



around T ~ 10 4 K for p > 10 -12 gem -3 . The prominent en- 
hancement of the opacity with increasing density is due to 
neutral hydrogen absorption, H + hv — > H + + e~ . 

The differences between the toy and full Pop III opac- 
ities are significant only at low temperatures. In the deep 
interior of the quasistar the opacity is constant, as it was 
for the analytic opacity, and we again find a convective 
zone that is described by an n = 3 polytrope for which 
p g /p r ~7.1mr ' 5 . 

In the outer region the flux is transported by radiative 
diffusion. The radiative layer covers between ~ 50 per cent 
and ~ 10 per cent of the quasistar by radius, becoming thin- 
ner as the photospheric temperature decreases toward T m i n 
(for a fixed jtibh) or as the black hole becomes more mas- 
sive (for a fixed T p h). In contrast, the percentage of mass in 
the outer radiative layer increases as the limiting temper- 
ature is approached. For ?71bh = 300 it goes from ~ 1 per 
cent at T p h = 10 4 K to ~ 18 per cent around T m i n , while 
for niBH = 7 x 10 3 it shows a similar behavior but never 
exceeds 2 per cent. 

The dramatic differences between structures computed 
using the real opacity, and those based on the toy opac- 
ity, are almost exclusively confined to the radiative layer. 
In Fig. [4] and Fig. [5] we show examples of the behavior 
of the opacity and the local Eddington ratio, Z/ ra d/£edd 
with Z/edd = 4nGM(r) c/n(r), for models computed with 
four different photospheric temperatures and constant black 
hole mass. The plots focus on the radiative zone. For clar- 
ity we plot temperature on the a>axis, since it spans more 
than an order of magnitude while the radius increases only 
by tenths of a per cent. As T p h decreases, the opacity in- 
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Figure 4. The opacity as a function of temperature in the outer 
layers of four quasistar models with Mbh = 300Mq and a = 0.1, 
computed usinf Pop III opacities. The photospheric temperatures 
Tph of the models are shown in the legends. On the right-hand 
panel we adopt a logarithmic scale for clarity, while the plot on 
the left-hand panel uses a linear scale. The temperature at the 
transition between the convective and radiative zones is denoted 
by the corresponding vertical lines. Note the rise in the bound-free 
peak at around T ~ 10 4 K (marked on the left panel by an arrow), 
which becomes dominant as the photospheric temperature drops 
and the photospheric density increases (see also Figs. [7] and Q5). 
Note also the steep drop in opacity at the photospheres. 



creases everywhere and the growth of the bound-free peak 
around T = 10 4 K becomes particularly prominent. Cor- 
respondingly, a peak in the I/ ra d/i e dd ratio forms. This 
peak first becomes super-Eddington for T ~ 10 K, and 
steadily grows and expands as the photospheric temperature 
drops until the entire radiative zone is super-Eddington at 
T P h = Tmin ~ 4 x 10 3 K. Note also that for T ph ~ 10 4 K a 
narrow super-Eddington region exists around T ~ 3 x 10 4 K. 
This is also true for higher photospheric temperatures. We 
comment on this point at the end of this section. 

Fig. U clearly shows the steep drop in opacity at the 
photosphere, and how it becomes more vertiginous as the 
quasistar cools down to T m i n . This feature ensures (as for 
the "toy" opacity) the presence of a minimum tempera- 
ture, below which no hydrostatic solutions can be found 
(Fig. [3J. As Tph decreases, the boundary condition at the 
photosphere (eq. [IS]) implies that the photospheric pressure 
passes through a minimum and then begins to increase. The 
photospheric radius and the transition radius R tI then stop 
increasing and start sinking rapidly inward in order to main- 
tain hydrostatic balance. Eventually no hydrostatic solution 
is possible. 

Despite the complicated behavior of the opacity, the 
analytic power-law scalings between quasistar mass, photo- 
spheric temperature, and black hole mass (eq. QT]) remain 
reliable at sufficiently high temperatures, although they suf- 
fer an offset in normalization correlated with the radial 
thickness of the radiative layer. To illustrate this, Fig. [5] 
shows the numerically computed m, vs T p h for models with 
Mbh = 300 Mq. The analytic scaling, 
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Figure 5. irad/^edd — 1 for the four models shown in Fig. [4] 
The solid horizontal line marks the border between the super- 
Eddington (above) and the sub-Eddington (below) regions. Note 
that for Tph ~ 10 4 K the super-Eddington zone is very narrow 
and confined around T ~ 3 X 10 4 K. For T ph < 10 4 K a peak 
appears around T ~ 10 4 K (marked by an arrow in the left-hand 
panel) and expands, eventually encompassing the whole radiative 
layer (dotted-line, right-hand panel). 
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Figure 6. The quasistar mass is plotted as a function of the pho- 
tospheric temperature (solid line) for models with mBH = 300 Mq 
and a = 0.1, computed with full Pop III opacities. The dotted 

20/7 
line shows the analytic scaling m, oc T . , obtained assuming 

that £tr = H = 1. The dashed line is also derived from the analytic 
relation given in eq. lilt , but with numerical values for £tr and 
k used instead of assuming that these factors are exactly unity. 
The star shows the mass of an (unphysical) numerical model com- 
puted with no radiative zone. This closely matches the analytic 
prediction, and indicates that the offset at high temperatures is 
due to the presence of a radiative zone. 



ocT 



20/7 



Ph 



(49) 
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tions (|36l) and ()39|) yields an expression for the density gra- 
dient, 



0.6 0.7 

r/R, 



Figure 7. Radial profiles of the interior temperature for the four 
models plotted in Fig. [4] For each model the transition radius Rtr 
is marked by the corresponding vertical line. 
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Figure 8. The same as Fig. [7] but for density. As T pn decreases, 
the radiative layer becomes denser. 



holds well down to T pn ~7x 10 3 K, though there is an off- 
set in the normalization. This offset is due to the presence 
of the radiative layer. A fully convective model at these high 
photospheric temperatures would lie on the analytic track: 
an example of such a model calculated for T p h = 9000 is 
plotted in Fig. [6] as a star. For higher marr, the offset is 
reduced as the radiative layer gets thinner. At lower tem- 
peratures, m* deviates slowly from the power-law until T m in 
is approached, whereupon to* drops more steeply. This is the 
catastrophic shrinking of the quasistar prior to evaporation. 
In this regime the analytic scalings do not hold. 

Although the integrated properties of the numerical 
Pop III opacity models resemble those derived analytically 
and numerically using the toy opacity, the presence of lo- 
cally super-Eddington fluxes for photospheric temperatures 
Tph > Tmin has substantial implications for the derived 
structures. Fig.[7]and Fig.[8]show the radial profiles of tem- 
perature and density for the models shown in Fig. [4] and 
Fig. [5] In hydrostatic equilibrium, the combination of equa- 



kT dp 
jjb dr 



GM{r) 



dp r 

dr 



pkdT 
/i dr 



(50) 



In regions where the radiative force —dp r /dr substantially 
exceeds the gravitational force the density gradient becomes 
positive and a local density inversion forms. At the same 
time, the temperature gradient steepens. In turn, the in- 
crease in the density and the decrease of temperature en- 
hance the contribution from bound-free hydrogen absorp- 
tion, raising the opacity further above n cs . This general 
behavior becomes more marked as the photospheric tem- 
perature drops, and the radiative layer becomes denser and 
cooler. 

The existence of density inversions in our models im- 
mediately raises the question of whether such structures are 
stable. There are two conceptually distinct concerns. In one 
dimension, the hydrostatic configurations we have computed 
are self-consistent, provided that the super-Eddington re- 
gions lie beneath the photosphere. An inward-directed force 
due to the gas pressure gradient compensates for the imbal- 
ance between radiation pressure and gravity. However, there 
could exist a second solution in which the super-Eddington 
flux drives a wind from the quasistar. If such solutions ex- 
ist, the hydrostatic solution could be unstable even in one 
dimension, leading to mass loss at temperatures well above 
T m i n . In two dimensions, it is even more doubtful that a 
density inversion would be stable. The resulting instabilities 
could lead to lateral density contrasts and non-magn etized 
photon bubbles of the type analyzed by [ShaviyJ (|200l|). 

We have not had notable success in analytically esti- 
mating the mass loss rate that might occur due to locally 
super-Eddington fluxes in the quasistar, and our numeri- 
cal scheme is not suited to tackle the problem. This fail- 
ure is hardly surprising, since the problem of what happens 
when a star develops a limited super-Eddington zone has 
been studied extensively in the context of Luminous Blue 
Variables such as n Carina (for a review, see Humphreys 
& Davidson 1994), apparently without any definitive theo- 
retical resolution being attained. One possibility is that a 
large-scale circulation pattern develops, superficially resem- 
bling convection, and the star suffers no mass loss at all 
(Owocki, Gayley & Shaviv 2004). Another possibility, em- 
pirically favored for LBVs, is that episodic mass loss occurs. 

For what follows, we conjecture that since the super- 
Eddington zones in our models arise due to high densities 
(and vanish if the density is reduced) , quasistars with Pop III 
opacities may develop outer regions in which there is circu- 
lation but little or no mass loss. As the width of the super- 
Eddington region becomes comparable to the radiative layer 
thickness, however, the extended acceleration zone is likely 
to permit a strong wind to develop. We therefore adopt T m i n 
as an estimate for the minimum temperature a quasistar can 
sustain before evaporating. That there are uncertainties in 
this identification should be very obvious. 
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4 CO-EVOLUTION OF BLACK HOLES AND 
QUASISTARS 

The properties of a quasistar will change as the black hole 
grows and, possibly, as the quasistar itself accretes matter 
from its environment. The thermal timescale in a quasistar 
interior is sufficiently short that the structure can adjust 
quasistatically. Therefore, we can model the co-evolution of 
the black hole and envelope as a series of equilibrium models, 
as long as the photospheric temperature exceeds the floor 
associated with the opacity crisis, T p h > Tmin. 

We will consider two models for the quasistar's inter- 
action with its environment. In the first, we will assume 
that the quasistar accretes matter at a constant rate, which 
we parameterize as O.lrho.iMo yr -1 . This growth rate is 
consistent with the BVR scenario for black hole growth in 
pregalactic halos. In the second, we will assume that the qua- 
sistar has a fixed mass, which we parameterize in units of 
10 6 Mq. We normalize time in years and assume, again mo- 
tivated by the BVR arguments, that the initial seed black 
hole mass is much smaller than the final mass attained prior 
to dissolution of the quasistar. 
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4.1 Before the Opacity Crisis 

As long as T p h > T m i n , the growth rate of the black hole is 
set by the Eddington limit of the quasistar, implying 



rhBH = 2.5 x 10 e .im*AfQ yr , 



(51) 



where e = O.leo.i is the accretion efficiency and we have 
taken l\ x = K — 1. The accretion efficiency relates the lu- 
minous output due to accretion with the rate of growth of 
the black hole mass. This efficiency factor is distinct from 
the Bondi efficiency factor a defined in eq. (0, and could 
be smaller than the standard value of 0.1 if energy is lost 
mechanically in a jet that does not couple to the enve- 
lope. When the quasistar grows at a steady rate, we have 
m» = O.lrho.ityr and the black hole mass grows according 
to 



m B H = 1.2 x 10 9 e 1 1 m .ity T 



1.2 x 10 _7 e^. 1 m^. 1 m*. (52) 



If the quasistar has a fixed mass, then the black hole grows 
linearly with time: 



niBH = 2.5 x 10 e .im*6*yr- 



(53) 



Equations (|11[) and (1521) imply that the photospheric tem- 
perature decreases as the black hole grows. The quasistar 
thus evolves toward the opacity crisis, reaching it when 
Tph = T m i n = 4000T m ,4 K. For the case of steady quasis- 
tar growth, this occurs when the quasistar mass reaches 



i o - in5 8/9 . -8/9 -4/9^,-20/9 

m*o = 1.8 x 10 e .i m o.i Q o.i T m,4 
and the black hole mass reaches 



^ 7/9 ■ -7/9 8/9 40/9 
fc 0.1 '"Ul "0.1 J m,4 



(54) 



(55) 



For a fixed envelope mass, the opacity crisis is reached when 
rriBH.o = 1.9 x 10 4 a^ /2 m^ 8 T- 5 / a . (56) 



Figure 9. Envelope mass versus black hole mass for static solu- 
tions at the minimum photospheric temperature. The dashed line 
is for the "toy" opacity model, the solid lines are for the numerical 
opacity. Static solutions are excluded in the lower shaded regions. 
The lighter shaded region is computed assuming that a = 0.1, the 
darker shaded region is for a = 0.05. Superimposed on the figure 
are the evolutionary tracks (eq. |52|). The upper track is for an 
accretion rate onto the envelope of M»= 1 Mq yr . The lower 
track is for M,= 0.1 M Q yr -1 . 



4.2 Numerical Results for Quasistar Evolution 

To test the above results, we numerically computed explicit 
evolutionary sequences by solving the equations, 



dAf BH 
dt 
AM, 



dt 



Tbh 
ec 2 



(57) 
(58) 



with Lbh(M*, Mbh,«) calculated from the model with no 
assumption as to the luminosity relative to the Eddington 
limit. The analytic evolution tracks for quasistars depend 
only on the assumption that the luminosity of the quasistar 
is close to the Eddington limit appropriate to the total mass. 
Our numerical results verified that this is a very good ap- 
proximation. Accordingly, to estimate the maximum mass a 
black hole can grow to within a quasistar it suffices to use 
the analytic growth track given by equation (|52[) . The ana- 
lytic estimate of the minimum temperature is less accurate. 
We therefore combine the analytic growth track with the 
numerically computed minimum photospheric temperature 
to derive the final mass. 

In Fig. [9] we plot analytic tracks corresponding to en- 
velope accretion rates of 0.1 Mq yr -1 and 1 Mq yr -1 . As 
the black hole grows, the tracks move toward the right-hand 
side of the plot, eventually crossing into the forbidden region 
set by the minimum photospheric temperature. The figure 
shows that the final black hole mass is higher for higher 
accretion rates onto the envelope, because the effective tem- 
perature decreases more slowly, delaying the final dissolution 
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at T m i n . We find that quasistars are indeed an efficient chan- 
nel for "growing" intermediate mass black holes. For a — O.f 
and rho.i ^ f , the final black hole mass is predicted to be at 
least a few thousand Solar masses. 



4.3 Scaling of Results with a 

The parameter a, which accounts for inefficiencies in the 
accretion flow within the Bondi radius, is largely unknown. 
The numerical results we have calculated are mostly for 
a — 0.1. Analytically, there is a simple scaling with a, valid 
for m, ^> 77iBH- In this limit, the only dependence of the 
quasistar structure on the black hole mass enters via eq. ([3]) 
for the luminosity, -Lbh oc Mb H q. Any solution in which the 
combination M — Mbh.Ch 1 ' 2 is constant should then have 
the same m» and T p h, along with the same radial structure 
well outside the Bondi radius. If a were to be smaller than 
our assumed value of 0.1, this would allow larger black holes 
to grow within the same envelopes. 

Numerically, the expectation that the results depend 
only on M is borne out for jtibh > 10 3 , irrespective of the 
photospheric temperature. For ttibh < 10 3 , it is valid only 
for photospheric temperatures well in excess of the mini- 
mum temperature. For low black hole masses and photo- 
spheric temperatures approaching the minimum value, the 
black hole mass cannot be ignored in the equation of hy- 
drostatic equilibrium, and the resulting structure depends 
separately on ttibh and a. In Fig. [9] we plot numerical re- 
sults for the minimum photospheric temperature computed 
with a = 0.1 and a = 0.05. As is evident from the figure, 
the offset between these curves is not constant. However, for 
the higher masses that are of interest when determining the 
maximum black hole mass that can be attained, the a 1 ' 2 
scaling is quite accurate. 



4.4 Post-Opacity Crisis: Dissolution of the 
Quasistar 

Our analysis of the junction between the radiative layer 
and convective zone indicates that no static solutions ex- 
ist for Tph < Tmin- In contrast to a red giant, protostar, 
or Thorne-Zytkow object, where feedback allows the energy 
source or envelope to adjust so that the photosphere fol- 
lows the Hayashi track, no stable feedback appears to ex- 
ist here. Once a quasistar reaches T m i n , the plummeting 
opacity causes the convective zone to release radiation at a 
super-Eddington rate. The deflation of the convective zone 
increases the rate of accretion onto the black hole, leading 
to a runaway. The only ways to avoid the destruction of the 
quasistar would be to decrease the mass of the black hole 
or to increase the mass of the envelope at an unrealistically 
high rate. 

A detailed analysis of the mass loss process is beyond 
the scope of this paper, but it is easy to show that the black 
hole is unlikely to grow much while the quasistar is evapo- 
rating. Once mass loss starts in earnest, the wind will carry 
away nearly all the energy released by the black hole. As- 
suming that the wind speed is of order the escape speed 
from the quasistar, vw ~ (GAf./fl,) 1 ' 2 , this implies that 
the quasistar evaporates when the energy liberated by ac- 
cretion equals the binding energy. (Although the quasistar 



is mainly radiation pressure-dominated, we assume that the 
binding energy is enhanced by rotation.) The mass accu- 
mulated by the black hole during the dissolution phase is 
then 



AM B 



GM 2 
eR* c 2 



on V 3 ™" 4 / 3 -7/8(77-4/3 B/f 
20e .i m o.i "o.oi T m,l m O, 



(59) 



which is negligibly small compared to 7TIbh,o. Therefore, we 
may regard mBH,o as the maximum mass attainable by a 
black hole growing inside a quasistar envelope. 



5 DISCUSSION AND CONCLUSIONS 

In this paper, we have studied the structure and evolution 
of quasistars, rapidly accreting black holes embedded within 
massive gas envelopes. We find that, for any black hole mass, 
there is a minimum photospheric temperature below which 
rapid dissolution of the envelope is inevitable. Both analytic 
and numerical models, computed using Pop III opacities, 
suggest that this minimum temperature is around 4000- 
5000 K. If quasistars are implicated in the formation of 
seeds for superm assive black holes in pre-gal a ctic h aloes, 
as suggested by iBegelman. Volonteri fc Reesl (12006) , this 
floor temperature implies that the most luminous quasistars 
would emit most strongly in the rest-frame near-IR. At typi- 
cal redshifts of z ~ 10, the observed spectrum would peak at 
A ~ 10/ttm. We defer discussion of the possible cosmological 
density of such sources, and hence their observability with 
future facilities such as the James Webb Space Telescope, to 
a subsequent paper. 

We have also studied the evolution of quasistars in sim- 
ple scenarios for the mass growth of their envelopes. Gener- 
ically, as the black hole mass grows the photospheric tem- 
perature falls, until eventually the limit imposed by the be- 
havior of the Pop III opacity is reached. Unveiling of the 
black hole appears inevitable long before it succeeds in ac- 
creting much of the envelope. Both analytic and numerical 
estimates suggest that seed black holes with masses between 
a few 10 3 Mq and 10 4 Mq are plausible outcomes of this sce- 
nario. The efficiency of the black hole accretion process that 
powers quasistars (parameterized here by a) is unknown — 
if the efficiency is low (due, for example, to polar outflows 
that couple poorly to the envelope) then larger seed black 
holes are possible. 

Although our numerical integrations of quasistar struc- 
ture confirm many aspects of the analytic model, they also 
reveal complex behavior in the radiative zone immediately 
beneath the photosphere. Narrow regions in which the flux 
is locally super-Eddington develop for temperatures signif- 
icantly in excess of the minimum temperature at which 
the entire radiative zone becomes super-Eddington. In our 
hydrostatic models, these zones are characterized by den- 
sity inversions. Interpreting this structure is tricky, since 
such density inversions may well be subject to one- or two- 
dimensional instabilities whose ultimate resolution is un- 
known. Hydrodynamic studies will be needed to determine 
the extent of mass loss that may occur at temperatures 
above the theoretical floor value. 

Finally, we note that here we have focused on black 
holes embedded within truly primordial gas at high redshift. 
Small amounts of metal pollution would act to increase the 
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opacity in the radiative zone, altering the minimum temper- 
ature and possibly increasing the likelihood of mass loss due 
to the formation of dust. Changes to the interior structure 
would be smaller. In particular, the density at the base of the 
radiative zone is so low (p ~ 10 -9 gem -3 ) that electron scat- 
tering will continue to dominate the opacity for T > 10 5 K, 
even in the presence of pollution. Of greater concern is the 
fact that metal enriched gas in the halo outside the quasistar 
would be more susceptible to fragmentation and star forma- 
tion. Too much star formation would reduce the rate of mass 
accretion onto the quasistar below the values (~ 0.1 Moyr -1 
or higher) that we have assumed, resulting in smaller final 
black hole masses. It is also possible that generically similar 
structures could form at lower redshift whenever small black 
holes encounter very high rates of gas inflow. Structures sim- 
ilar to those we have described could allow stellar remnants 
at the center of merging galaxies to grow significantly via a 
transient quasistar stage. 
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